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SAMPLED-DATA TECHNIQUES APPLIED TO A DIGITAL 
CONTROLLER FOR AN ALTITUDE AUTOPILOT 
By Stanley F. Schmidt and Eleanor V. Harper 


SUMMARY 


Sampled-data theory, using the Z transformation, is applied to the 
design of a digital controller for an aircraft -altitude autopilot. Par- 
ticular attention is focused on the sensitivity of the design to parameter 
variations and the abruptness of the response, that is, the normal accel- 
eration required to carry out a transient maneuver. Consideration of 
these two characteristics of the system has shown that the finite settling 
time design method produces an unacceptable system, primarily because of 
the high sensitivity of the response to parameter variations, although 
abruptness can be controlled by increasing the sampling period. Also 
demonstrated is the importance of having well-damped poles or zeros if 
cancellation is attempted in the design methods. 

A different method of smoothing the response and obtaining a design 
which is not excessively sensitive is proposed, and examples are carried 
through to demonstrate the validity of the procedure. This method is 
based on design concepts of continuous systems , and it is shown that if 
no pole-zero cancellations are allowed in the design, one can obtain a 
response which is not too abrupt, is relatively insensitive to parameter 
variations, and is not sensitive to practical limits on control-surface 
rate. This particular design also has the simplest possible pulse 
transfer function for the digital controller. 

Simulation techniques and root loci are used for the verification of 
the design philosophy. 


INTRODUCTION 


The design of sampled-data systems, the theory of which is applicable 
to feedback control systems utilizing a digital computer, is a relatively 
new field but is receiving considerable attention in those applications 
in which the basic information is received as pulses or in the form of 
numbers. This may be due to greater accuracy requirements, since digital 
transducers can be made more accurate than analog transducers, or to the 
fact that the actual measurement device, for example a radar, yields a 
sampled signal. Another application in which digital computers in real 
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time control systems are receiving attention is where a large amount of 
flexibility of changing computation procedures is required (e.g., an 
interceptor aircraft having a fire control system, an automatic landing 
system, altitude and direction hold systems, etc.). In these applica- 
tions the digital computer offers a possibility of weight reduction and 
improved reliability, both highly desirable features. 

The sampled-data theory 1 has progressed to the point where it is 
possible to design a digital-controller pulse transfer function so as to 
obtain an over-all system pulse transfer function which will meet certain 
specifications, such as being stable and physically realizable, having a 
minimum settling time in response to a given input, having a ripple -free 
response after a short transient, etc. A survey of the literature, how- 
ever, has indicated a scarcity of information on the sensitivity of a 
system to parameter variations or how to control, in the design process, 
the abruptness of the response. References 2 and 3 touch upon the sub- 
ject of abruptness and show one method of smoothing the response by 
proper design of the characteristic equation. Another shortcoming in the 
literature has been that the examples are restricted to plants of third 
order or less. 

It is the purpose of this paper to show, first, that considerations 
of the sensitivity of the response to parameter variations and abruptness 
of the response to transient inputs preclude the use of finite settling 
time design (i.e., one in which the error is reduced to zero in a speci- 
fied number of sampling instants); second, :o show that continuous system 
design concepts can be used to select the dominant poles of the closed- 
loop pulse transfer function so as to achieve a desired transient 
response; and third, to present results of a study of the design for a 
fourth-order plant. 

In order to obtain some indication of vhe practical limitations 
imposed on the response of the system, the ’plant" is taken to be an 
aircraft and the system to be one that controls the aircraft's altitude. 
The abruptness and sensitivity of the response to a transient maneuver 
will be judged in this application by the magnitude of the maximum load 
factor and the change in the stability of the system with parameter 
variations . 

In order to obtain a system response that will be acceptable to a 
variety of inputs and flight conditions the theory of design used here 
is different from that normally employed to "soften" the response. 


^he theory and design methods contained in this report were derived 
to a large extent from a course taught at Stanford University by G. F. 
Franklin. Reference 1 by J. R. Ragazzini and G. F. Franklin contains an 
extensive bibliography of the available literature on sampled-data theory. 
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NOTATION 


SLry 




normal acceleration, ft/sec 2 
mean aerodynamic chord, ft 

pitching-moment coefficient, 
vertical -force coefficient, 


pitching moment 
(l/2)pV 2 Sc 

vertical force 
(l/2)pV 2 Sc 


2.7183 


h altitude, ft 

Ah change in altitude, ft 

Iy moment of inertia, lh-ft/sec 2 


j 

\/ "I 

m 

mass, lb sec 2 /ft 

Ma, Mq 

1 pV* 2 Sc /5C m ^C'm c 

2 Iy \ ba ’ bqc/2V 2V/ 

% 

1 pV 2 Sc / dC m jr ^C m \ 

2 Iy \dcxG / 2 V 2V ' 35/ 


pitching velocity, radians/; 

s 

Laplace operator 

S 

wing area, ft 2 

T 

sampling period, sec 

V 

velocity, ft/sec 

z 

s T 

e = Z transform operator 


1 pVS fic z dCz\ 

Z a> Z 8 2 m ^ ^ ' 36/ 


Z[f(t)] Z transform of f(t) 
a angle of attack, radians 
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7 

5 

i 

e 

n 

P 

u>n 

F(z _1 ) 


flight path angle , radians 
control-surface deflection , radians 
damping ratio 
pitch angle , radians 


3 . 1^159 

density , lb sec 2 /ft 4 
natural frequency , radians/sec 


n 


finite polynomial in z 1 of forn ^ a-j_z 1 


1=0 


Superscript 


* sampled version of a function of time 


FUNDAMENTAL CONCE :5 TS 


The theory of sampled-data systems is covered in considerable detail 
in references 1 and 4 . However, it is desirable to include herein certain 
of the fundamental concepts of sampled-data systems so that the applica- 
tion of the theory to the example chosen wi_l be understandable by those 
familiar with the design and analysis of linear feedback systems. 


Elements of a Sampled-Da:a System 


The basic elements of a simple sampled-data system are shown in 
sketch (a). Here the continuous error signal e(t) is measured at 
regular intervals, T seconds apart, by the sampler. The output of the 



Sketch (a) 
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sampler consists of a train of pulses whose amplitude (or area) represents 
the value of the input at each sampling instant. This pulse train e*(t) 
is modified by the controller to provide stability and other characteris- 
tics to the complete system as well as smoothing of the sampled data. 

The controlled- system characteristics are presumed known and in this 
report will be represented by a transfer function relating altitude to 
controller output. 


Analytical Description of the Sampler 

The output of the sampler shown in sketch (a) can be written as 

e(t)i(t) = e*(t) (1) 

where i(t) is a train of unit impulses occurring every T seconds. 

This operation is illustrated in sketch (b) 



Sketch (b) 


Now using an infinite series expression for i(t) the sampler output 
becomes 

00 

e*(t) = e(t) ^u Q (t - nT) (2) 

n=o 

where u 0 is the unit impulse function. Since the unit impulse function 
is zero except at time nT, this equation can be rewritten as 

00 

e*(t) = Ve(nT)u 0 (t - nT) (3) 

n=o 

where e(nT) is the value of the input when t = nT. The Laplace trans- 
form of this equation is 



6 


E*(s) = ^e(nT)e n ^ s 
n=o 

A second method of specifying E*(s) is by expanding 

+ 00 

^Tu c (t - nT) = i(t) 

- 00 

in a Fourier series which gives 

+00 

j 2 rtnt/T 
e 

-OO 



OO 


(5) 


since the Fourier series coefficients are constant and all equal to l/T; 
therefore , 


By theorem (see ref. 5) 


then, 


e*(t) = 

e(t) 

+ 00 

1 V j 2 Jtn 1 ./T 

Tl 6 
- 00 

(6) 

5) 






L 

f(t) 

jut - 

e 

= 

F(s - t u>) 

(7) 



+ 00 



E 

*(s) 

-*> 

El 

G - ¥) 

(8) 


This formula shows that E*(s) is a periodic function repeating itself 
every 2 tc/T radians per second as illustrated in sketch (c). 



T T 


Sketch (c) 
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The Z Transform 


The Z transform of a function of "time (defined only ah sampling 
instants) can be found by substituting z = e s ^* in equation ( 4 ). Thus 

00 

Z[e*(t)] = E(z) = ^e(nT)z' n (9) 

n=o 


The infinite summation can be found in closed form for all applica- 
tions considered here. To find the closed form one can either sum the 
series ( 9 ) or use the complex convolution integral 


E (z ) = E*(s) 


n c + joo 


z = e 


sT 




E(A) 


C - joo 


1 -e 


dA 

-T(s - A) 


( 10 ) 


|z = e 


sT 


Tables have been prepared for a large number of the useful Z transforms 
(refs. 1 and 4 ); thus, in general it is not necessary to use complex 
integration. The Z transform can relate the sampled input of a system 
to the sampled output. When it is used in this manner it is commonly 
referred to as a "pulse" transfer function. 


To obtain e(nT) from a closed 
that can be taken 


1 . 

Refer to tables 

(refs . 

2. 

Expand 

E ( z ) in 

powers 

3 - 

Use the 

formula 

e(nT) 


the unit circle. 


Z transform there are three courses 
1 and 4 ) . 

of z 1 by long division. 

= 2jtJ f E ( z ) z31 1(iz where T is 


The initial and final value of a pulse sequence which results from 
inversion of a pulse transfer function can be determined readily from 
the following theorems: 

Initial value theorem f(o) = lim F(z) (n) 

Z — ► 00 

Final value theorem f(oo) = lim (l- z _1 )F(z) (12) 

z-*i 

A further important point about the Z transform is the stability of a 
given pulse transfer function. Since z = e sT it can be seen that the 
left half of the s plane maps into a unit circle in the z plane. 

Thus , if all poles of a pulse transfer function are located inside this 
unit circle the pulse transfer function is stable. 
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In general, for a sampled-data system of the type illustrated in 
sketch (a), that is a sampling of the error signal, the Z transform can 
he used similarly to the Laplace transform for a continuous system. 
Root-locus methods of analysis can be used in the z plane; however, 
they are difficult to apply since the region of interest is that inside 
the unit circle. The drawing of the root locus becomes relatively tedious 
since asymptotic behavior provides no help in giving a rough idea of the 
position of the loci as is the case in the s plane. It should also be 
recognized that the Z transform design specifies only what happens at 
sampling instants. This is not a serious drawback, however, for in order 
to determine system behavior it is always possible to test the designed 
system by simulation techniques. Also, from the transfer function of the 
continuous elements of the system, a relatively accurate idea of the 
system behavior between sampling instants can be determined. 

Because of the drawbacks of root -locus techniques in the z plane 
and for other reasons, the design of sampled-data systems is probably 
best carried out by a different method. This method, which is explained 
in the next section, could be applied to continuous system design as well; 
however, to the authors' knowledge, it is not In widespread use. 


Design Criteria For Digital-Controller Pulse 
Transfer Function 


The procedures for the design of the digital controller D(z) are 
outlined with reference to the simplified blcck diagram shown in sketch (d). 
In this sketch D(z) represents the pulse transfer function of the digital 
controller. The zero-order hold circuit constructs a continuous signal 



Sketch (d) 
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from a sampled signal as shown in sketch (e). The transfer function for 
this operation is (l-e" sT )/s. The block labeled "plant" represents the 
fixed elements of the system. 




( 13 ) 

( 14 ) 

( 15 ) 

(16) 


Equation (l6) is the fundamental design equation. Certain mathematical 
constraints must be put on this equation in order to arrive at a stable 
system which has the desired characteristics. These constraints and the 
reasons for them are as follows : 


I. Stability 

A. K(z) must be of the form 
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[ (l - a 1 z“ 1 )(l - a 2 z x ). . • ](p c , + P 1 z -1 + . . .) ( 17 ) 

where the £ f s are undetermined coefficients and z = a. ± , 
z = a 2 , . . . are the zeros of G(z) on or outside the unit 

circle (on or outside the stable region in the z plane). 

The reason for this constraint is obvious from equation (l6) 
If it is not satisfied D(z) will have a pole where G(z) 
has a zero, and consequently some of the poles of D(z) will 
be outside the stable region causing the over-all system to 
be unstable since perfect canc ellation can never be expected 

B. l-K(z) must be of the form 

[(1 -biz'hO- -baz' 1 ). • .](7o + 7 1 z" 1 + . • • ) (18) 

where the 7 *s are undetermined coefficients and z = b x , 
z = b 2 , . . . are the poles of G(z) on or outside the unit 
circle in the z plane. The reason for this constraint is 
also obvious from equation (16), for if it were not met, 
zeros of D(z) would be imperfectly canceling poles of G(z) 
outside the stable region, a fact which would cause insta- 
bility in the over -al l , system. 

II. Zero Steady -State Error to an Input of the Form r(t) = t n 

The Z transform of t n is of the form [f(z -1 )]/[ (l - z" 1 ) 11 * 1 ] . 
If zero error to such an input is desired (at sampling instants) 
then 1 - K(z) must be of the form 

(1 - z" 1 ) n+1 [F 1 (z’^ U] (19) 

where Fi(z _1 ) is a polynomial in z” 1 satisfying other 
constraints. This can be seen by applying the final value 
theorem (12) to equation (l4). 

III. Transient Performance 

A. Finite settling time 

(1) Minimum 

If only I and II are satisfied then K(z) will be the 
lowest order polynomial in z -1 . This results in a 
minimum finite settling tune since reference to equa- 
tion (9) shows that e(nT) will be zero (to a unit 
pulse input r(t)) after T times the order of K(z) 
seconds . 

(2) Zero ripple 

Zero ripple by definition means that for an input of 
form r(t) = t n the output must be of the form t n . 
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In other words, the output must follow the input between 
sampling instants. With reference to sketch (d) this 
can only be accomplished if e^nT) is zero or a constant 
in the steady state, for otherwise the output of the 
zero-order hold circuit would be stepping from one value 
to another, which would cause ripple. 


Ei(z) 


C(z) 

G(z) 


R(z) 


K(z) 

G(z) 


( 20 ) 


For e 1 (nT) equal to a constant in the steady state, 
equation (20 ) must be of the form 

F(z~ x ) 

1 - z" 1 


as can be seen from equation (12). Thus, for zero 
ripple all the zeros of G(z) must be contained in K(z). 
This is not the only requirement. In addition if 


R(z) = 


F(z~ x ) 

_ -1 \ n +l 


then G(z) must have an nth-order pole at z = 1 to 
cancel all but one of the poles of R(z). An even 
simpler way of assuring that zero ripple is possible is 
to count the number of poles of G(s) at s = 0. If 
r(t) = t n then the number of poles of G(s) at s = 0 
must be greater than or equal to n + 1. 


B. Smoothing 


As is described in reference 1 it is possible to add a 
denominator 2 to K(z) which has been determined by the pre- 
ceding results so that K(z) is of the form 


K(z ) = 


F(z~ 1 ) 

(l + c 1 z“ 1 + c 2 z -2 + 


The addition of this denominator changes the system such 
that finite settling time is no longer obtained. It is 
mentioned in reference 1 that the choices of c’s can be 
made for 


^ote that in reality we are not necessarily increasing the order of 
the characteristic equation. What is really done is to change the charac- 
teristic equation from one with all roots at z = 0 to one which allows 
some of the roots to be at other positions inside the unit circle in the 
z plane . 
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(1) Smoothing the transient response so that it is not 
as abrupt. 

(2) Optimizing the response in the presence of noise. 

This report illustrates that if considerations of sensitivity 
and abruptness of response are important, then one should 
always design the system so that K(z) has a denominator. 

It is shown in a later section of this report that a choice 
of c’s in equation ( 21 ) which will satisfy abruptness of 
response requirements can be determined from continuous 
system design theorems. 


DESCRIPTION OF PROBLEM 


The problem chosen as an illustrative example is an altitude command- 
type autopilot. In order to reduce the problem to block diagram form it 
is necessary to derive the aerodynamic equations and put them into trans- 
fer function form. In order to do this certain assumptions must be made 
as follows: 

(a) The aircraft does not roll. 

(b) The velocity is a constant. 

(c) The altitude is approximately constant. 

(d) Small angle approximations are valid for 6 and 7 . 

(e) The aircraft is initially in trim flight. 

(f) The aerodynamic coefficients are constant. 

With these assumptions the following equations are valid 



-q+d = -7 = ZqCX + Z^S 

(22) 


q = MqO t + Mqq + Mg 5 *r M^d 

(23) 

q = 0 

pitching velocity 


N 

sJ 

II 

> 

normal acceleration 


II 

•N- 

> 

rate of change of altitude 



from which the transfer functions given below can be derived: 
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1. (Mg + M^Zg)s+ (M^Zg MgZ,-^) 

5 " + ^ + 

Yl = a* = ( -Z S )s 2 + (Mg + M^Zgs - (MgZ a - F^Zg ) 

5 8 s2- (Za. + Mq+MaJs + CMqZa-Ma) 


(24) 

(25) 


As a representative aircraft an interceptor at a cruising velocity of 
695 feet per second, 30,000 feet altitude was chosen. The aerodynamic 
constants are tabulated below 

Mg = -l6. 04 
Mg, = -O.163O 
z a = -0.6716 
V = +695 ft/sec 

With these values the transfer functions (24) and ( 25 ) become 

q. _ -l 6 . 02 s - 9.961 

5 _ s 2 + 1.324s + 7.065 


Mq = -0.489 
% = -6.737 
Zg = - 0.1205 


vj a z 83 . 75 s 2 + 54 . 6 ls - 6923 
s ” 6 s 2 + 1.324s + 7.065 


(27) 


The block diagram of the system is as shown in figure 1. Two cases 
designed for comparative purposes are as follows: 

Case I Kq - 0 

Case II Kq = -0.18 

These values of Kq result in the transfer functions between M, the 
output of the hold circuit, and a z of 


Case I 


5 

M 


= 1.0 


fz f^z 

M 5 


fz = 83 . 75 s 2 + ^4,6is - 692 3 
M s 2 + 4.206s + 8.858 


(28) 


Case II 
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APPLICATION OF DESIGN CRITERIA TO ALTITUDE AUTOPILOT 


Theoretically, the only constraint which must be imposed, upon the 
design of a digital controller is that the system be stable. As shown 
in the previous section a minimum finite settling time system would 
result from this constraint alone. Since it is desirable to have zero 
error continuously one must also impose the zero ripple constraint. 
Practical considerations, however, such as abruptness of the response to 
step inputs, smoothing of noisy inputs, and sensitivity to parameter 
changes may require that additional constraints on the design be imposed. 
It is desirable, therefore, to investigate the different designs which 
can be made and study their characteristics by means of simulation. 


Choice of Sampling Period 


There seems to be no theoretical method of choosing the optimum 
sampling period for a closed-loop system. If the input command were 
band limited to a frequency u>, then the sampling theorem states that if 
l/2u> = T, all information can be reconstructed from the samples. This 
theorem does not apply here since the system is completely satisfactory 
only if zero error to all possible inputs is maintained. This dilemma 
is not studied here. The choice of sampling period is based on knowledge 
of the control of an aircraft by a human pilot. Previous studies indi- 
cated that a period of 0.25 second is about ihe longest that should be 
chosen. For exemplary purposes a period of 3 second is also chosen in 
order that effects of sampling period in the design can be demonstrated. 
Practical considerations probably would force the choice of period to be 
the shorter of the two, principally because £. 1 second sampling period 
would cause a rough flight for the pilot. The 0.25 period or 4 cps 
sampling rate appears to be high enough that the jerkiness of the hold 
circuit output during transient inputs would be quite well filtered by 
the control-surface servo and aerodynamic lags. 


Finite Settling Time Design With Zero Ripple 


Finite settling time design, Case I, Kq = 0, T = 0.25 second.- The 

steps involved in the design are shown for this example for the benefit 
of those unf ami liar with digital -controller c.esign methods. The first 
step is the determination of G(z) 


G(s) 



/83.75s 2 + 54.6 :.s - 69235 . /j\ 
V s2 + 1.324s +7.065 ) V S V 


( 29 ) 


Transfer function ( 29 ) is expanded in partial, fractions 
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G(s) 



979-8 191-3 114.7 

S 3 + S 2 + S 


114,7s + 343.1 \ 
s2 +1. 324s + 7.0657 U ' 


The individual terms are converted to Z transforms, using the tables 
of reference 1 or 4, and then recombined to give 


G(z) 


1.360z 3 -12.67z 2 -11.96Z + 1.074 
(z - 1) 2 ( Z 2- I.356Z +0.7182) 


(3D 


= 1.360z 1 - 12.67z~ 2 - 11.96z~ 3 + 1.074z~ 4 , 

(l - z _1 ) 2 (l - l.356z _1 + o.7l82z“ 2 ) 

The next steps are to determine K(z) and l-K(z) such that stability 
and zero ripple constraints are satisfied. 


1 


” K(z ) 


K(z ) = (Numerator of G(z) in powers of z _1 )(p o + p i z -1 ) 

= (1.36oz- 1 -12.67z- 2 -h.96z- 3 + i.074z- 4 )(p o + 3 1 z- 1 ) (33) 

Poles of G(z) on or outside unit]/ -1 -2 -3, 

= L circle in powers of z^ j^o+r^ +7 2 z +7 3 ^ ) 


= (1 -2z _1 + z" 2 )(7 0 +7 1 z“ 1 + 7 2 z" 2 +7 3 z" 3 ) (34) 

The coefficients of equal powers of z _1 in the expression l-K(z) as 
obtained from equations (33) and (34) are equated. Enough undetermined 
coefficients must be provided to give a sufficient number of simultaneous 
equations in y's and p f s to allow their solution. For this example, the 
solution of the simultaneous equations gives 


7 0 = 1.000 
7 X = 2.215 
y 2 = 1.276 
7 3 = -0.1211 

Then 

K(z) = (l.360z -1 -12. 67z -2 -11. 96z“ 3 + 1.074z -4 )(-0. 1579 + 0. H28 z _1 ) 

(35) 

1 -K(z) = (1 - z _1 ) 2 (l + 2. 215z _1 + 1.276z" 2 -0.12Hz -3 ) (36) 

The last step is the substitution of equations (35), (36 ), and (32) into 
equation (l6) to give 


Po = -0.1579 
P x = 0.1128 



1 6 


-0.1579(1 ~ O.Tl^z" 1 )^ - l.3^6z~ 1 + 0 . 7 l 82 z~ g ) 
1+ 2.215z _1 + 1 . 276 z ' 2 - 0.1211z" 3 


The design at this stage is complete. This system as designed, 
however, may he unsatisfactory because of sensitivity to parameter 
changes and abruptness of the response. Ar analog computer simulation 
is one of the best ways of analyzing such a system. This method is used 
here. Appendix A shows one procedure by which a digital computer can be 
simulated on an analog computer. This method is derived in reference 6. 
Derived, also in appendix A, is the detailed computer diagram for this 
example. Appendix B contains a description of the electronic sample- 
hold circuit used for simulation purposes in this investigation. 

The results of the simulation studies are summarized in figure 2. 

There are two important things to note in figure 2. The first is that 
the system does not have a finite settling time. This is obvious from 
figures 2(a) and (b) in which it can be seen that both the step response 
and the ramp response contain an oscillatoiy mode which damps exponen- 
tially. This appears to be due to the extieme sensitivity of the system 
which makes it impossible to simulate the system even with highly accu- 
rate analog computing equipment. The secord item of importance to note 
is that for a step input of 10 feet, figure 2(a), a peak acceleration of 
approximately 31 g T s is required during the transient maneuver. This 
large peak in acceleration is a result of fhe linear analysis. In an 
actual aircraft, control-surface position end rate limits as well as 
aerodynamic nonlinearities would prevent such an excessive peak. One 
cannot state what the exact performance of the finite settling time 
design would be with the actual nonlinearii ies; however, the introduction 
of a high-performance servo in the simulation resulted in an instability 
of the system which indicates that in all probability nonlinearities 
would also cause system instability. 

Figure 2(c) illustrates the effects of a 10-percent increase and 
decrease in system gain. These results plus the fact that the simulated 
system does not have finite settling time indicate that the design is 
quite sensitive to parameter changes. 

Figure 3 is a root -locus plot showing effects of gain on theoretical 
pole locations of a closed-loop system. Tie system is noted to be very 
sensitive to gain changes around the desigred gain which places the closed- 
loop poles at the origin. 

Figures 4(a) and (b) are root-locus plots showing effects of altitude 
on the closed-loop pole-zero locations. Ncte that a drop in altitude 
changes all the aerodynamic coefficients because of the consequent change 
in air density. Again it is seen that the pole locations which are at 
the origin for 30,000 feet move a considerable distance (two almost 
becoming unstable) for 27,000 feet. 
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In s umma ry, then, for two reasons - (a) relatively high sensitivity 
to parameter changes and (b) response which is much too abrupt - the 
finite settling design method does not result in a satisfactory closed- 
loop performance for this particular case. 

Finite settling time design. Case II , Kq = -0.18, T = 0.25 second. - 

The results of the previous system studied showed that the theoretical 
finite settling time design could not be simulated. Figures 3 and 
illustrate the relatively high sensitivity of the system; however, this 
is not necessarily the only reason that the simulated system did not 
agree with theory. Another reason might be related to the fact that 
zeros of D(z) theoretically must cancel poles of G(z). Because the two 
canceled poles have a relatively low damping ratio (£ = 0.258), it is 
possible that imperfect cancellation, in addition to the high sensitivity 
previously noted, could be an important factor in the disagreement. The 
artificial damping added in Case II increases the damping ratio of the 
canceled poles to 0-707; thus making it possible to study the effect of 
damping on the ability to simulate a finite settling time design. 

For this example, the use of the same procedures as those of the 
previous section gives the following pulse transfer functions : 


G(z) 


1.019z 3 ~ 10.21z 2 - 7 . 527 z + 0.7295 
(z -l) 2 (z 2 - 1.022z + 0.3 J *-9 1 0 


(38) 


K(z) = (1.019a -1 -10.21z -2 - 7.527z" 3 + 0.7295z" 4 )(-0.2154 + 0.1528z“ 1 ) 

(39) 

l-K(z) = (1 - z _1 ) 2 (l +2.220z _1 + 1.085z" 2 -0.1115Z -3 ) (4o) 

_ -0.215M1 ~ 0-7096z~ 1 )(l -1.022 z~ 1 + 0.3494z~ a ) 

D(Z ' l + 2.220z _1 + 1 . 085 z ‘ 2 - 0.1115z“ 3 

The results of the simulated response for this example are summarized 
in figure 5- It can be noted from both figures 5(a) and 5(b) that the 
simulated system has a finite settling time, a fact which agrees with the 
theory. Figure 5(a) shows this system also to be very abrupt, as would 
be expected, requiring about 33 g's peak for a 10-foot step input. Fig- 
ure 5(c) illustrates that the response of this system is also quite 
sensitive to gain variations . 

From the results of the simulated systems. Case I and Case II, it 
can be concluded that if the design method being utilized requires pole- 
zero cancellation one should be certain the canceled poles of G(s) are 
in a well-damped region of the s plane. 
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Finite settling time design, T = 1.0 second . - The two cases studied 
previously are unacceptable for application to an aircraft due to both 
sensitivity to parameter variations and abruptness of response. One of 
the questions it is necessary to answer for a sampled-data system is: 

Do these two parameters vary together or can they be independently con- 
trolled? It should be obvious that one way of controlling abruptness of 
response is by increasing the sampling period. Thus, one can maintain 
a finite settling time design and reduce the abruptness of response, but 
the important question is will the sensitivity to parameter c han ges also 
be reduced? 

An alternative procedure which could be utilized is to keep the 
sampling period the same but specify, as another constraint, the time 
response at sampling instants to be smooth and as slow as desired for as 
many sampling instants as is necessary. If one maintains a finite 
settling time design in this manner the result will be that the order of 
the system is increased by one for each sampling instant specified; thus, 
the digital computer will become more and more complicated. This second 
method will not be illustrated here because of its additio nal complexity. 

To determine the effect of sampling period on abruptness and 
sensitivity, the two cases are designed for finite settling t ime period 
T = 1 second. For Case I, Kq = 0, the pulse transfer functions are: 

G(z) - 7^.9s -TB.H .. . 

(z -1) 2 (z2 + 0.8704; + 0.2662) ^ ' 

K( z ) = (-162. 8z -1 - 1104 z ' 2 - 754. 9z~ 3 - 72.11 5 -4 ) (-0.001601 + 0 . 001123 Z" 1 ) 

(^3) 

i-k(z) = (l -2z _1 + z -2 )(i +i.739z _1 -f 0 . 8947 Z " 2 +o.o8io6z" 3 ) (44) 


, ) = -0.001601(1 -0.7016z~ 1 )(1+0.8704z~ 1 +0.2662z~ 2 ) 


This system was simulated by the method of appendix A and typical 
transient responses are summarized in figure 6. Again it may be noted 
that in the simulation, finite settling time was not achieved. The 
abruptness of the response has been considerably reduced, that is, peak 
acceleration is about O.87 g’s for a 10-foor step. Figure 6(c) shows 
the response to be somewhat sensitive to ga:.n changes; however, the best 
indication of a sensitivity problem is the inability to simulate a theo- 
retically finite settling time design. The reason for this could be the 
attempt to cancel poles of G(z) which are in relatively low damped 
regions as was the case for T = 0.25 secona. This particular system is 
too oscillatory to be useful so we shall nov consider Case II where 
canceled poles are more heavily damped. In this case the pulse 
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transfer functions are: 


-nooz' 1 - 572.lz~ g - 200. 6z~ 3 - 7.468z~ 4 
(1 - z- 1 ) 2 (l + 0.1247z~ 1 + 0.0l491z- 2 ) 


( 46 ) 


K(z) = (-110.3Z -1 -572.1Z -2 - 200. 6z“ 3 - 7-468z -4 ) (-0.003501 + 0.002378z _1 ) 

(47) 

l-K(z) = (l - z -1 ) 2 (i + i.6l4z~ 1 + o.4866z~ 2 +0.01780z -3 ) (48) 

Df , = -0.003501(1 -0.6793 z~ 1 )(i + 0.1247z- 1 + o.oi491z~ 2 ) 

l + l.6l4z _1 + o.4866z~ 2 +0.01780Z -3 ^ 

The results of this simulation are summarized in figure 7* Fig- 
ures 7(a) and (b) show that finite settling time was achieved for this 
case, and illustrates again that if poles of G(z) are to be canceled 
by zeros of D(z) then these poles must be heavily damped. Figure 7(a) 
also illustrates, as was anticipated, that the abruptness of response 
has been reduced considerably being about I .30 g's for a 10-foot step 
input. Figure 7(c) also gives conclusive evidence that the finite 
settling time design is quite sensitive to gain variations even though 
the abruptness of response has been reduced. Further evidence of this 
sensitivity problem is illustrated in figure 8 which is a root-locus plot 
of Case II for open-loop gain variations. 

The general summary of the results would tend to confirm the fact 
that a finite settling time design procedure will always result in a 
sensitive system and, thus, if it is to be used for any practical appli- 
cations, the plant must have a transfer function free of nonlinearities 
and with coefficients that are nonvariant. An intuitive reason for this 
is illustrated by the root-locus plot of figure 8. We see that a finite 
settling time design requires all poles of the closed-loop system, with 
the exception of those being canceled, to be at the origin in the z 

It can be seen the origin is a position of extreme sensitivity 
since open-loop gain changes of ±5 percent cause the poles to move a 
considerable distance; whereas, an additional change of ±5 percent causes 
only the small motion indicated in the figure. The alternative procedure 
for reducing the abruptness of the response, previously mentioned, would 
place an even greater number of poles at the origin. From the results 
shown in figure 8, one would believe that placing more poles at the origin 
would increase the sensitivity rather than decrease it. 

It must be concluded from these standpoints that some of the poles 
of the closed-loop system should be located in other positions than the 
origin inside the unit circle of the z plane if we are to obtain satis- 
factory performance. The next section illustrates a method for choosing 
this location based on continuous -system analogy. 
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Considerations Involved in Adding a Denominator to K(z) 


The two previous systems studied indicate that it would "be desirable 
to reduce the sensitivity to parameter variations and smooth the tran- 
sient response. It is mentioned in references 2 and 3 that the addition 
of a denominator (sometimes referred to as a staleness factor if the 
denominator is first order) will smooth the transient response. It will 
be shown here that if the denominator is properly chosen, the sensitivity 
to parameter changes is reduced. 

There are two possible methods for selecting a denominator. The 
first is a simple trial and error process which can be very tedious. 

The second is to base the sampled-data system response on what might be 
reasonable for a continuous system of the s;une type. By means of the 
latter method, experience gained from similar continuous system designs 
can be used to select the proper location of the dominant modes in the 
s plane which are known to give satisfactory performance. These modes 
or location of poles can then be transferred to the z plane and the 
design carried through using the added denominator. 

Consider the continuous system shown ii the block diagram of fig- 
ure 9. The problem is to design the network D(s) so that the over-all 
closed-loop response h Q /hi will be satisfactory. Figure 10 shows the 
pole and zero locations of G(s) for the tw:> cases. It is obvious that 
for Case I the complex poles of the aircraft are insufficiently damped. 

For this case, D(s) must either (l) have complex zeros to attract the 
aircraft poles to a more favorable position or ( 2 ) cancel the poles with 
zeros and place new poles in a more favorable position. The second choice 
can never be used for an aircraft since (l) poles and zeros are never 
known very accurately, (2) their positions shift with Mach number and 
altitude, and ( 3 ) gust inputs would excite the oscillatory mode since it 
is not canceled for inputs other than those from the control surface. 

Figure 10(b) illustrates that the effect of adding an inner-loop 
pitch-rate feedback (fig. l) is to shift the poles to a more favorable 
position. It should be noted that aircraft automatic control systems 
almost always utilize an inner-loop autopilct of either the combined 
normal acceleration and pitch-rate feedback type or the simple pitch-rate 
feedback illustrated in this example. The addition of normal accelera- 
tion as a feedback along with pitch rate alLows both the natural frequency 
and damping of the aircraft modes to be shifted substantially from the 
basic airframe oscillatory mode. 

It should be noted that D(s) for the continuous system (Case II ) 
could be a simple lead network provided the zero, pole, and gain of the 
network were chosen so that the complex poles did not shift to an 
unfavorable position for the closed-loop performance. Case II is a much 
more satisfactory system for the continuous system since the network can 
be simpler. It will therefore be used for the following sampled-data 
system studies . 
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For Case II then, what might be a reasonable selection of the domi- 
nant second-order mode of the continuous system? This question cannot 
be answered without specifying the task the system is required to perform. 
Since this report is dealing principally with the application of sampled- 
data design techniques, the only concern will be that the closed-loop 
system be relatively insensitive to parameter changes. The dominant mode 
of the continuous system will be the closed-loop location of the two poles 
which are at s = 0 open loop. Figure 11 illustrates a root-locus plot 
for one location of the zero and pole of a simple lead network for D(s). 
Note that the zeros of G(s) are omitted since their effects on the loci 
near the origin are negligible , except that the gain must be negative 
for stability. It should be noted that there is not much change in the 
aircraft oscillatory mode if the gain is varied from 0 to 17. A nominal 
operating gain for the system was chosen as 17 which corresponds to the 
dominant mode at u>n = 0.5 radian per second and £ = 0.707. From the 
figure it can be seen that the resultant over-all closed-loop transfer 
function can be approximated by 

hp _ 4s + 1 

h i _ 4s 2 + 2 . 828 s +1 


(50) 


This approximation is justified because the other three poles and two 
zeros are a relatively long distance from the dominant real zero and 
complex-pole locations. It will be assumed here that this transfer func- 
tion, given by equation (50), satisfactorily performs the task for which 
the continuous system is being designed. 

The dominant mode (denominator) of the sampled-data system closed- 
loop, K(z), will thus be chosen at the z plane location of the two 
poles of the denominator of the transfer function given by equation ( 50 ). 
This transforms (for T = 0.25) into 

(z -0.9126+0. 8l7i)(z -O. 9126 - 0- 08171) = z 2 -1.825 z + 0.8395 (51) 

Figure 12 illustrates the location of poles and zeros of G(z) for Case II, 
Kq = -0.l8, T = 0.25. It should be noted that the two open-loop poles at 
z = 1 will move along the dotted lines as the gain is increased to arrive 
at the desired location given by equation (51). It should also be noted 
that D(z) can be chosen for this case in two different ways. The first 
method, cancellation permitted, is to cancel the two complex poles of the 
aircraft by zeros and place other poles so that, at the desired gain, they 
end up at the origin. The second method, cancellation not permitted, is 
to choose a single pole for D(z) and force K(z) to have at least four 
poles away from the origin. This second method will actually show that 
the digital-controller pulse transfer function D(z) can be simplified by 
a properly chosen denominator. Both of these methods will be studied in 
the following sections. 
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K(z) 


[Fi(z" 1 ) > meeting other constraints] (p Q + (3 1 z” 1 + . . .) 
Desired denominator 


(52) 


and 


[F 2 (z“ 1 )^ meeting other constraints] (y Q + + . . . ) 

Desired denominator 

(53) 

By subtracting equation (52) from unity and equating coefficients of like 
powers of z -1 to those of equation (53) "the P's and y's are uniquely 
determined. For this example all zeros of G(z) appear in K(z) and the 
two poles at z = 1 appear in l-K(z) 


G(z) 


1.019Z- 1 -10.21z~ 2 -7.52Tz~ 3 + 0.7295z~ 4 
(1 - z _1 ) 2 (l - 1.022z _1 + 0.3494z -2 ) 


( 54 ) 


K(z) 


( 1 . 019 Z" 1 - 10.21z~ g - 7.527z ~ 3 + 0.7g95z~ 4 )(p o d^z" 1 ) 
(1 - l.825z _1 + 0 . 8395 z“ 2 ) 


(55) 


1 -K(z) = h - z ~ 1 ) 2 (/o + 7 iZ~ 1 + 7 2 z ~ 2 + 7 3 z~ 3 ) 

(1 -1.825z" 1 + 0.8395 z -2 ) 

The solution for the y ' s and p 1 s gives 

y 0 = 1.000 

7 1 = 0.18730 

7 2 = 0.07768 

7 3 = -0.0082657 

Use of equation (l6) gives 

. . _ - 0 . 01223(1 - 0 . 9264 z~ 1 )(l -1.022z~ 1 + 0 . 3494 z~ 2 ) 
z “ 1 + 0 .l 873 z _1 + 0 . 077672 " 2 -0.008266Z" 3 


Po = -0.01223 

Pi = 0.01133 


(56) 


(57) 


The closed-loop sampled-data system was simulated by the same method as 
that of previous studies. The results of the simulation are summarized 
in figure 13 . As can be noted in figure 13(a), a 10-foot step now calls 
for a little over a 1 g maneuver. The response is much slower than for 
that of the finite settling time design. The abruptness of the control- 
surface motion at the beginning of the transient causes this peak in 
normal acceleration and can be traced in part to the design method of 
canceling the poles of the aircraft. Figure 13(c) illustrates that the 
system response is not very sensitive to gain variations. 
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The effects of adding a simulated control- surface servo of transfer 
function 


5 _ 1 

M 1 + 0.02s 


(58) 


is shown in figure 13 (d) for three values of control-surface rate limit- 
ing, y . The inclusion of a linear servo alone gave identical response 
to figure 13 (a) so that the large overshoots are directly traceable to 
effects of control- surface limiting. It is interesting to note that the 
overshoot is as high as 500 percent for §max = 20° per second. The 
unusual fact that the overshoot is higher for Smax = 20 ° per second than 
for &max = 10 ° per second is due to the time relationship between the 
actual control-surface motion and the command input during the initial 
part of the transient. The reason for the overshoot is traceable to the 
design method of canceling poles with zeros of D(z). This might be 
expected since any nonlinearities such as saturation in the control sur- 
face servo velocity make cancellation impossible for a transient input. 


Cancellation permitted. Case II, Kq = -0.l8, T : 1.0 second . - The 
application of the method previously described for this example gives 


-110. 3z -1 - 572. lz ~ 2 - 200 . 6 z ~ 3 - 7.468z~ 4 
(1 - z _1 ) 2 (l +0.1247z _1 + 0.0l491z“ 2 ) 


(59) 


K(z) = 


(- 110 . 3 z -1 - 572 . iz ~ 2 - 200 . 6 z ~ 3 - 7 . 468z~ 4 )( -0.0009842 + 0.0007864z~ 1 ) 


1 - K(z ) = 


1 -1.320Z -1 + O.4966Z -2 


(1 


.- 1^2 


(60) 


) g (l + 0 . 5710 z " 1 + 0.l622z ~ 2 +0.00587z~ 3 ) 

l - l . 320 z" 1 + 0 . 4966z“ 2 ' 


-0.0009842(1 - 0 . 7990 z _1 )(l +0.1247z" 1 + 0.0l491z- 2 ) . 

D(z) = ; 7 5 o (62) 

1 + 0.5710Z 1 + 0.l622z 2 + 0.00587z 3 


Simulation of this system by the method described in appendix A gave the 
results which are summarized in figure lk. Of particular note in compar- 
ing figures 13 (a) and l4(a) is that increasing the sampling period does 
not particularly increase the response time even though there is a con- 
siderable reduction in both control -surface deflection and maximum normal 
acceleration. Figure l4(c) illustrates that the response is not particu- 
larly sensitive to parameter changes. The effect of 6 limiting on this 
system was studied; however, for a 10 -foot step, the called for control- 
surface motion is so small that control-surface rate limiting produced 
no noticeable effect. 
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Cancellation not permitted. Case II, Kq = -0.l8, T = 0.25 second. - 

As was previously mentioned the digital -controller pulse transfer func- 
tion D(z) can be simplified by choosing a different characteristic 
equation. It may be noted in figure 12 that only a single zero and 
single pole of D(z), if their location is properly chosen, will result 
in a stable design. One of the questions, then, is how to locate the 
pole and zero. Root-locus methods were tried; however, for a s ample d- 
data system which has as many poles and zeros as this case the root-locus 
method is very tedious. 

An alternative method was used as follovs: 

K(z) = [Fi(z- 1 )](P 0 + 3 1 z~ : -+. • .) 

' (Desired dominant mode ) (l + c ; z -1 + c 2 z“ 2 + . . .) ' 


l-K(z) 


[F 2 (z X )] (7q + 7 ! s ‘ ~ 1 + - . . ) 

(Desired dominant mode ) (l + cqz" 1 + c 2 z -2 + . . .) 


(64) 


where F 1 (z _1 ) contains all the zeros of G(s ) and F 2 (z -1 ) contains all 
the poles of G(z). For this example, Case ]I, Kq = -0.l8, T = 0.25 


K(z) 


(1.019Z- 1 -10.2lz~ g - 7.527z~ 3 + 0.7295z- 4 )(p o + PiZ" 1 ) 
(1 - 1.825Z -1 +0.8395z _2 )(l + C 1 z~ 1 + c 2 z -2 ) 


(65) 


(1 - Z~ 1 ) 2 (l -1.022z~ 1 + 0.3494z~ 2 )(y o +7 1 z~ 1 ) 
(1 -1.825z-i + 0.8395z“ 2 )(1 + c 1 z-i + c 2 z-2) 


(66) 


Solution of simultaneous equations for p ! s, /'s, and c*s gives 


P 0 = -0.002763 7o = 1-000 
P x = 0.002588 = -0.005402 

By use of equation (l6) then 


c x = -1.205 
c 2 = 0.4018 


D(z) = 1 0-002763(1 - 0.9365Z-U 

1 - 0.005402z _1 v j 

A root -locus plot for this design is illustrated in figure 15 . It 
may be noted that the closed-loop complex poles (due to the aerodynamics) 
move to a somewhat more highly damped position in the z plane. Their 
motion is not large, however, and it would be reasonable to assume that 
the transient performance will be relatively insensitive to position 
variations. This conclusion is demonstrated in the simulated responses 
summarized in figure l6. 
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Figure l6(d) illustrates that this design which allows no cancella- 
tion is insensitive to reasonable values of control-surface limiting. 

The differences in output are not even large enough to be recorded. The 
only noticeable difference is in the initial control-surface motion which 
i is completed within one sampling instant and, for the input magnitudes 

tested, created no noticeable effect on the output. 

Figure 16(e) shows that effects of a change in the aerodynamics 
corresponding to a change in altitudes of 10,000 and 20,000 feet are 
relatively negligible. For this case the D(z) was left unchanged at 
the design value of 30,000 feet altitude. This insensitivity is princi- 
pally due to the relatively "low-gain" system. Had a more rapid response 
been required the same conclusions might not have been reached. 

It is of interest to compare the transient responses of the sampled- 
data system derived by this design procedure with those of the continuous 
system which was designed to obtain desired dominant mode characteristics. 
The root -locus plot of the continuous system used for comparison is shown 
in figure 11. Figures 17(a) and (b) show the very close comparison of the 
transient responses of the two systems. It should be noted that the start 
of the sampled-data system response was shifted so that the two start at 
the same point since a delay of 0 to 0.25 second (dependent on the time 
of the application of the step) can be experienced by the sampled system. 


CONCLUDING REMARKS 


It has been demonstrated that the finite settling time design of a 
sampled-data system does not produce desirable characteristics when the 
method is applied to the design of an altitude autopilot. Two undesirable 
features are the very abrupt response and the extreme sensitivity to 
parameter changes which are unsatisfactory for aircraft and many other 
automatic control systems. Increasing the sampling period generally 
reduces the abruptness of response; however, it does not appear to elimi- 
nate the sensitivity problem. The sampling period may also be dictated 
by other considerations, so that it might not be a variable. 

It has been demonstrated that basing the sampled-data system dominant 
mode on a continuous system design results in a much more practical system 
with respect to the two previous considerations. It ban also been shown 
that cancellation of aircraft poles by their inverse in the stabilizing 
digital network is undesirable. This, of course, is not surprising since 
the same conclusion is true in continuous systems. 

The most satisfactory design method tested was one in which no 
cancellation was allowed and the resulting system has the simplest digital 
control function possible. 
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In general, the use of the Z transform for sampled-data system 
design appears to be almost as easy to handle as the Laplace transform 
for continuous system design. One drawback which has been noticed is 
that root-locus methods are not easily applicable, principally because 
the Z transform of the open-loop system has almost the same number of 
zeros as poles. Since stability requires poles to be confined to the 
unit circle, a knowledge of the asymptotic behavior for a large z does 
not help. In addition a solution for saddle points in the z plane is 
alm ost impossible. This means that one must use the tedious method of 
trying points to see if they are on the loci without knowing the approxi- 
mate loci. The other drawback is that for the high order system of this 
example the sampling period should be chosen initially rather than carried 
through as an arbitrary constant to be deterndned in the final design. 

This is simply because of the complexity of the resultant equations which 
occur if it is arbitrary. Since there does rot appear to be a clear cut 
means of selecting the sampling period, then one would, in general, have 
to duplicate several designs at different periods in order to arrive at 
the most suitable one for the particular application. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Jan. l4, 1999 
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APPENDIX A 

SIMULATION OF THE SAMPLED-DATA SYSTEM 


A sampled-data system can be simulated on a conventional analog 
computer. The only additional device required is a sample-hold circuit. 
There are many ways of constructing this type device utilizing relays or 
various electronic circuits. The circuitry used for this study is briefl; 
described in appendix B and is an all electronic device. 

Given a sample-hold circuit, how is D(z) simulated? The method 
used is described in reference 6; however, the technique is simple and 
will be described here. The following equations can be developed from 
figure 18 


A(z) = 


A(z)B(z) . 

D(z) = = A(z)B(z ) 

Z[(l-e- sT )/s] 

Z[(l-e- sT )/s] 


1+ Z {[ (1 -e~ sT )/s]Q(s)} 1+ (l - z -1 )Z[Q(s)/s] 

from equation (A2) 


~Q(s)l 


r 1 

i 

s J 


La(z) x j 

d-z-n 


B(z) = Z 

f r om equat ion ( Ak ) 


1 - e 


-sT 


P(s) 


1 - e 


-s 


= (l-z^JZ 




'P(s)' 


B(z) 

(1 - z-i) 


(Al) 

(A2) 


(A3) 


(A4) 


(A5) 


A given D(z) is to be simulated by the circuit of figure 18. The prob- 
lem is then to split D(z) into A(z) and B(z) so that Q(s) and P(s) 
can be determined from equations (3) and (5)> respectively. 


D(z) = K 


(l -a 1 z~ 1 )(l - a 2 z~ g )(. 
.(1 -baz-i)(l -b 2 z-2)(. 



(A6) 
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Now the following conditions insure that Q(s) and P(s) will be physically 
realizable stable networks: 

(1) All unstable poles of D(z) shall be contained in A(z). 

(2) Unstable zeros of D(z) shall be contained in B(z). 

(3) All gain, K, shall be contained in B(z). 

(4) Stable poles and zeros shall be distributed between A(z) 
and B(z) to allow { [ l/ A ( z ) ] - l}[l/(l-z“" 1 )] and B(z)/(l-z“ 1 ) 
to be expanded in partial fractions of type (3/ (l - az"' 1 ) . 

Reference 6 gives a number of examples and explains in detail how all 
physically realizable D(z) f s may be constricted from R-C networks. 

Since the desire here is to simulate, one may allow complex poles. Con- 
ditions 1, 2, and 4 assure Q(s) and P(s) will be stable networks. These 
conditions certainly do not appear to be necessary for simulation pur- 
poses, a fact which essentially means that all of D(z) except the gain, 

K, can be assigned to A(z). In figure l 8 , this can be seen to be a 
desirable way of handling the simulation since only one sample-hold 
circuit is required. For the example, Case I, Kq = 0, T = 0.25 


D(z) = -0.1579 


(l - 0.7l46z -1 )(l - 1 . 356 z _1 + o. 7182 z- 2 )' 

1 + 2.215Z" 1 + 1.2f6z -2 -0.1211z -3 


A(z) 


B(z) = -0.1579 = Gain of D(z) 

(1 -0.7l46z~ 1 )(l- 1.356 s" 1 + o. 7 l 82 z -g ) 

1 + 2.215z _1 + 1.276z" 2 -0.12Hz -3 


(at) 

(a 8 ) 

(A9) 


Q(s) 


La(z) 


-1 


l 


4.285Z -1 - 0.4-lllz -2 + 0.394lz -3 
(l-z -1 ) “(1-Z -1 )(1 - 0.7l+6z -1 )(l -1.356 z -1 + 0.7182z -2 ) 

(A10) 


This expression is expanded in partial fraction expansion 


Z 


Q(s) 

s 


= z 


41.26 30.83 


io.4oz + 0.8049 


- 1 z -0.7146 “ z 2 -1.356z + 0.7182 


(All) 


The terms are individually converted to the s plane and recombined to 
give 


_ 8.655s 2 +84.385 +392.0 

“ (s +I.345)(s2 + 1. 324s + 7 . 065 ) 


(A12) 
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It should be noted that Q(s) has a lower order numerator than denomina- 
tor. This is absolutely essential because of the circuit used and is a 
consequence of imposing condition 3 . 

A block diagram of the system to be simulated is shown in figure 19 . 
Figure 20 shows the corresponding analog-computer diagram used for 
simulating the system. 



30 


APPENDIX B 

AN ELECTRONIC SAMPLE -HOLD CIRCUIT 


As was previously mentioned, there are a number of methods of 
constructing a sample-hold circuit. The circuitry utilized for this 
example is all electronic and for some applications may offer advantages 
over circuitry utilizing relays. 

Figure 21 is a block diagram of the circuit used. The diode bridge 
acts as a gate circuit connecting the output of amplifier 1 to the input 
of integrator 2. When a positive pulse is applied to A and a negative 
pulse to B, the diode bridge connects the output of amplifier 1 to the 
input of integrator 2. When zero signal is applied to A and B the input 
to integrator 2 is open; thus, the diode bridge acts as a sampler of 
amplifier 1. The function of integrator 2 is to hold the output between 
the sampling instants. Integrator 2 acts as an open-loop integrator dur- 
ing the time in which the pulse is off and, as a consequence, its output 
(other than for drift) stays at the value it was at the previous sampling 
instant. If an initial step of Ej_ is assumed the circuit can be seen 
to work in the following manner. Initially amplifier 1 builds up to the 
step value; as soon as the pulse arrives at A and B, the integrator 
output starts to change value. If the time constant of amplifier 1, 2, 

3 combination is about l/5 of the sampling time, then the output E 0 
will be equal to +Ej_ at the end of the sample time. Then E 0 is held 
at Ei by the integrator condenser. Drifts are corrected each sampling 
instant. For this problem it was desired than the sample time should not 
be longer than 1 percent of the sample period, T, and that the problem 
be run on a real time basis on the analog computer. This means that for 
a 4 cps sampling frequency the sample time * h/lOO X 0.25 second = 2.5 
milliseconds and the time constant of amplifier 1, 2, 3 combination 
should be approximately 0-5 millisecond. 

As can be seen from the above number, relatively fast amplifiers 
must be available for use in this sample-hold circuit. This is one 
difficulty which would probably, in general, jriake this circuit difficult 
to fabricate from conventional analog compute- amplifiers. For use in 
this problem specially designed high-perf ormaice chopper-stabilized 
amplifiers were available. No particular problems were noticed in the 
use of this circuit other than the need for ai occasional adjustment to 
compensate for drift in the integrator. 
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Figure 2.- Transient responses for a finite settling time design; 

Case I, T = 0.25. 
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(a) Pole location. 

Figure 4 .- Effect of altitude on closed loop pole and zero location 
Case I, T = 0.25, finite settling time design. 
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(b) Zero location. 
Figure 4 .- Concluded. 
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Figure 5 ■ " Continued 
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.ock diagram of a continuous altitude -command autopilot 
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Figure 10.- Location of poles and zeros of G(s). 
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continuous system; Case II. 
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(c) Effects of ±10-percent gain changes; 10-foot step input. 

Figure 13 .- Continued. 
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(d) Effects of a rate-limited, control-surface servo; 10-foot step input 


Figure 13 •- Concluded. 
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(a) 10-foot step input. 


Figure 17- - Comparison of transient responses of a sampled-data system 
cancellation not permitted, with a continuous system. 









Figure 19.- Block diagram for a finite settling time design; Case I, 

T = 0.25. 
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Figure 20.- Analog computer circuit for a finite settling time design; Case I, T = 0.25. 




71 



Figure 21.- Block diagram of an electronic sample hold circuit. 
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